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1. Introduction 

In a series of papers, we and others have re- 
ported the resuhs of global fits to pre-LHC [iHS] 
and LHC 2010 data [7HTU] in the frameworks of 
simplified variants of the minimal supersymmet- 
ric extension of the Standard Model (MSSM) [11] 
with universal supersymmetry-breaking mass pa- 
rameters at the GUT scale. We consider a class 
of models in which R-parity is conserved and the 
lightest supersymmetric particle (LSP), assumed 
to be the lightest neutralino Xi [H], provides 
the cosmological cold dark matter [13]. The spe- 
cific models studied have included the constrained 
MSSM (CMSSM) [H], with parameters mo, mi/2 
and Aq denoting common scalar, fermionic and 
trilinear soft supersymmetry-breaking parame- 
ters at the GUT scale, and tan/3 denoting the 
ratio of the two vacuum expectation values of 
the two Higgs fields. Other models studied in- 
clude a model in which common supersymmetry- 
breaking contributions to the Higgs masses are 
allowed to be non-universal (the NUHMl) jl5j . 
a very constrained model in which trilinear and 
bilinear soft supersymmetry-breaking parameters 
are related (the VCMSSM) [TB] , and minimal su- 
pergravity (mSUGRA) [I6lfl7] in which the grav- 
itino mass is required to be the same as the uni- 
versal soft supersymmetry-breaking scalar mass 
before renormalization. 

The impressive increase in the accumulated 
LHC luminosity combined with the rapid analyses 
of LHC data by the ATLAS [TSlfTO]. CMS [20H22] 
and LHCb Collaborations [25 is putting increas- 
ing pressure on these and other supersymmetric 
models, in the continuing absence of any signal for 
supersymmetry. In this paper we update our pre- 
vious frequentist fits [8l[9] to include the analyses 
of ^ 1 /fb of LHC data made public in July and 
August 2011 at the EPS HEP and Lepton-Photon 
Conferences, termed here LHCi/fb, and also dis- 
cuss the impact of the result on BR(i?s M^/^ ) 
by the CDF Collaboration [24]. As in [8], we also 
incorporate the results of the direct search for 
dark matter scattering by the XENONIOO Col- 
laboration [25] 0. 

^Preliminary versions of these results were posted on [9] 
on July 25th, 2011. 



The approach we use has been documented in 
our previous papers [THS1[71[5], so we do not de- 
scribe it in detail here, concentrating on rele- 
vant new aspects. We construct a global like- 
lihood function that receives contributions from 
the standard portfolio of electroweak precision 
observables, as well as B-decay measurements 
such as BR(6 — s- 57) and BR(i?„ — ;> ri/r)- 
The contributions to the likelihood function from 
BR{Bs ->■ M+M"), the XENONIOO direct search 
for dark matter scattering and the LHC searches 
for supersymmetric signals are calculated within 
the MasterCode framework [S]. Concerning the 
theoretical predictions for the different observ- 
ables, the MasterCode framework incorporates a 
code for the electroweak observables based on [26] 
as weU as the SoftSUSY FeynHiggs [28], 

SuFla [29], Superlso [30], MicrOMEGAs [3l] and 
SSARD [32] codes, using the SUSY Les Houches 
Accord [33] ■ The MasterCode framework is such 
that new observables can easily be incorporated 
via new 'afterburners', as we discuss below for the 
LHCi/fb constraints. We use a Markov Chain 
Monte Carlo (MCMC) approach to sample the 
parameter spaces of supersymmetric models, and 
the results of this paper are based on a basic re- 
sampling of the CMSSM with TOM points and 
a resampling of the NUHMl with TOM addi- 
tional points, both extending up to ?tt.0j™i/2 = 
4000 GeV. 

Our results are based on the public results of 
searches for supersymmetric signals using ~ 1/fb 
of LHC data analyzed by the ATLAS and CMS 
Collaborations and ^ 0.3/fb of data analyzed 
by the LHCb Collaboration. For our purposes, 
some of the most important constraints are pro- 
vided by the ATLAS [H] and CMS [20] searches 
for jets -|- l^T events without leptons, as well 
as searches for the heavier MSSM Higgs bosons, 
H/A [T9II2T]. Also important are the new upper 
limits on BR{Bs M+M") from the CMS [12], 
LHCb [ig and CDF Collaborations [H], which 
we incorporate in this paper as described be- 
low H. In this paper we focus on the analysis 
of the effects in the CMSSM and the NUHMl. 



^For other studies of recent data on BR(_Bs — > fJ,'^ fJ- ) in 
supersymmetric frameworks, see 1341 . 
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As discussed briefly below, the VCMSSM and 
mSUGRA models are further disfavoured by the 
LHCi/fb data. 

The absences of signals in the jets + I^t 
searches disfavour the ranges of the model mass 
parameters {mQ,mi/2) that had been favoured 
in our previous analyses of the CMSSM and 
NUHMl [HIS], and our current best flts have 
Too 150 to 450 GeV and mi/2 750 GeV. 
Reconciling these larger values of (toq, rni/2) with 
{g — 2)f^ favours values of tan/3 ^ 40, though with 
a large uncertainty. The regions of parameter 
space with large tan /? are constrained also by the 
new upper limits on BR(i?s ^'^fi~), as well as 
the LHC H/A searches. Using our standard im- 
plementation of the ((7 — 2)^ constraint based on a 
Standard Model (SM) calculation [42j, and com- 
bining the theoretical and experimental errors in 
BR(6 S7) in quadrature, we find that the p- 
value for the CMSSM best-fit point is now ~ 15%, 
and that for the NUHMl is ^ 16%. On the other 
hand, we show that if the {g — 2)^ constraint is 
dropped much larger regions of the (too,toi/2) 
and other parameter planes are allowed at the 
68 and 95% CL, and these p-values increase to 
38% in both models. In contrast, changing the 
treatment of BR(6 sj) by adding linearly the 
theoretical and experimental errors has relatively 
little impact on the fits and their p-values, in- 
creasing them both to 18%. 

On the basis of these results, we present 
updated predictions for the gluino mass mg, 
BR{Bs — > and the fight and heavy Higgs 

masses Mh and Ma- We also present updated 
predictions for the spin-independent dark matter 
scattering cross section, cTp', stressing the impor- 
tance of the uncertainty in the 7r-nucleon a term 
S^AT [35]. 

In addition, we use our results to present like- 
lihood functions for the thresholds for sparticle 
pair production in e+e^ collisions. These results 
indicate that, within the CMSSM and NUHMl, 
the best-fit values for the sparticle thresholds lie 
above Ecu = 500 GeV. However, we empha- 
size that these results are derived in the con- 
text of specific models with specific universal 
soft supersymmetry-breaking masses at the GUT 
scale, and do not apply to other classes of super- 



symmetric models. 

2. Implementations of the New LHC Con- 
straints 

Jets +]^T searches 

The CMS and ATLAS Collaborations have 
both announced new exclusions in the (toq, mi/2) 
plane of the CMSSM based on searches for events 
with jets -|- J^T unaccompanied by charged lep- 
tons, assuming tan/3 = 10, Aq = and /i > 0. 
The updated CMS ax analysis is based on 1.1 /fb 
of data [20], and the updated ATLAS 0-lepton 
analysis is based on 1.04/fb of data [18]. It is 
known that 0-lepton analyses are in general rela- 
tively insensitive to the tan /3 and ^0 parameters 
of the CMSSM, as has been confirmed specifically 
for the CMS ax analysis, and they are also insen- 
sitive to the amount of Higgs non-universality in 
the NUHMl. Therefore, we treat these analy- 
ses as constraints in the (TOo,mi/2) planes of the 
CMSSM and NUHMl that are independent of the 
other model parameters. The ATLAS [36] and 
CMS Collaborations [37] have also announced 
new exclusions for searches for jets -I- I^t events 
with one or more charged leptons with ~ 1/fb 
of data, but these have in general less expected 
sensitivity, and are more dependent on the other 
model parameters, so we do not include them in 
our analysis. A similar remark applies to the new 
ATLAS limits on events with b jets -|- I^t unac- 
companied by charged leptons using 0.83/fb of 
data [38] and on events with 6 jets + + 1 
lepton using 1.03/fb of data [39] B 

The CMS and ATLAS 0-lepton searches are 
more powerful in complementary regions of the 
(mo,mi/2) plane. Along each ray in this plane, 
we compare the expected CMS and ATLAS sen- 
sitivities, select the search that has the stronger 
expected 95% CL limit ^and apply the constraint 
imposed by that searcliQ We assign Ax^ = 5.99, 

■^It would facilitate the modelling of LHC constraints on 
SUSY if each Collaboration could combine the results from 
its different missing-energy searches, as is already done for 
Higgs searches. 

^It would also facilitate the modelling of LHC constraints 
on supcrsymmetry if the results from different Collabora- 
tions were combined officially, as was done at LEP, is al- 
ready done for BR(i3i, — > fj,^fi~) searches, and is planned 
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corresponding to 1.96 effective standard devia- 
tions, along the CMS and ATLAS 95% 0-lcpton 
exclusion contours in the plane. In 

the absence of more complete experimental in- 
formation, we approximate the impact of these 
constraints by assuming that event numbers scale 
along rays in this plane oc Ai^"^ where A4 = 
^TOq + rn^^^ (more details can be found in ^). 

We then use these numbers to calculate the ef- 
fective numbers of standard deviations and cor- 
responding values of Ax^ at each point in the 
plane. 

Searches for heavy MSSM Higgs bosons 

The CMS Collaboration has announced a new 
constraint on the heavy MSSM Higgs bosons from 
a search for the neutral bosons H/A — > r+r" us- 
ing 1.6/fb of data [21] and ATLAS has presented 
a similar constraint using 1.06/fb of data [12], the 
results being presented as 95% CL upper limits on 
the product of the production cross section times 
T+r" branching ratio as a function of the com- 
mon mass, Mh/Aj for masses smaller than about 
500 GeV. In our analysis we use the CMS con- 
straint, which has the greater expected sensitiv- 
ity. The CMS Collaboration has also announced 
a constraint on the decay chain t — !> iJ+ — r+z^ 
using 1.09/fb of data [401, but this yields a con- 
straint in a generic {Ma-, tan /3) plane that is much 
weaker than the above searches for H/A, so we do 
not implement it in our analysis. 

We assign Ax^ = 3.84, corresponding to 1.96 
effective standard deviations, to model parame- 
ter sets predicting an H/A signal at the 95% CL 
given by the CMS constraint, for each fixed value 
of Ma- Other model parameter sets are assigned 
values of Ax^ according to the numbers of ef- 
fective standard deviations corresponding to the 
numbers of events they predict. For any fixed 
value of Mjj/ji, these event numbers scale approx- 
imately as (tan/3)^. 

Constraints on BR(_Bs — > ^i^) 

Three new results on BR(_Bs — > jj,^ ji^) have 
been announced recently. One is an excess of 
candidate Bs ^ events reported by the 



for Higgs searches. 



CDF Collaboration [53], which corresponds to 
BR(B, ^ = (1.8lJ-j) X 10-8 or ^R{B, ^ 

fi+l-i ) < 4.0 X 10"*^ at the 95% CL. The other two 
new results are upper limits from the CMS Col- 
laboration using 1.14/fb of data [22]: BR{Bs 

< 1.9 X 10-8 at the 95% CL, and from the 
LHCb Collaboration using 0.34/pb of data [23]: 
BR{Bs M'^M") < 1-5 X 10-8 at the 95% CL. 

These three results are reasonably compatible, 
though there is some tension between the CDF 
and CMS/LHCb results. The two latter coUab- 
orations have released an official combination of 
their results [H], which yields BK{Bs fi^ 
< 1.08 X 10-8 at the 95% CL. In our implemen- 
tation of the BR(i3s — > constraint, we use 
Ax^ corresponding to the full likelihood function 
provided by this combination, which has a global 
minimum close to the SM prediction. We also 
comment on the changes in our results that would 
follow from an (unofficial) combination with the 
CDF result [24], which would yield a Ax^ func- 
tion with a minimum at BR(_Bs — at 
about twice the SM value. 



Constraints on dark matter scattering 

We incorporate the upper limit on the spin- 
independent dark matter scattering cross sec- 
tion erf provided by the XENONIOO CoUabo- 
ration |25j in a similar manner to [8]. In that 
paper we discussed extensively the uncertainty in 
the spin- independent scattering matrix element 
induced by the relatively ill-determined value of 
the TT-nucleon a term, S^rAf- In this paper we use 
S^AT = 50±14 MeV, and neglect other uncertain- 
ties, e.g., in modelling the dark matter distribu- 
tion. We also do not consider here other experi- 
ments reporting signatures that would require (7^^ 
above the XENONIOO limit. Nor do we consider 
limits on spin-dependent dark matter scattering 
and astrophysical signatures of dark matter an- 
nihilations, which currently do not constrain the 
CMSSM and NUHMl [5]. 
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3. Results 

The {mo, 1111/2) planes in the CMSSM and 
NUHMl 

Fig. [T] displays contours with Ax^ ~ 2.30 (red) 
and 5.99 (blue) relative to the minimum values 
of at the best-fit points (denoted by green 
stars) in the (mo, mi/2) planes for the CMSSM 
and NUHMl Such contours are commonly in- 
terpreted as 68 and 95% CL contours. The solid 
lines are the contours after incorporation of the 
LHCi/fb results, and the dotted lines are the CL 
contours obtained from an analysis of the pre- 
LHC and pre-XENONlOO data 0. Consequently, 
the differences in the contours show the full im- 
pact of the ^ 1/fb data set of LHC data. The 
crinkles in these contours give an indication of 
the sampling uncertainties in our analysis. 

We see that the new best-fit points with 
(too, mi/2) = (450, 780) GeV in the CMSSM and 
(150, 730) GeV in the NUHMl (denoted by solid 
green stars) lie well within the previous 95% CL 
region. On the other hand, the pre-LHC best- 
fit points with (too,TOi/2) = (90, 360) GeV in 
the CMSSM and (110,340) GeV in the NUHMl 
(denoted by open stars), lie far outside the re- 
gions allowed by the LHCi/fb data. Thus, we 
see that there is now significant tension between 
the LHCi/fb and pre-LHC data sets. The full set 
of parameters of the post-LHCi/fb and pre-LHC 
best-fit points are shown in Table [T]|f|. 

In both the CMSSM and the NUHMl, as we 
discuss later, the 68% CL upper limits on TO2/2 in 
Fig. [Dare largely driven by {g — 2)^, and the 95% 
CL upper limits are largely driven by the relic 
density constraint, with large-TOi/2 points lying 
in the heavy-Higgs rapid-annihilation funnel at 
large tan/3. The fact that much larger regions in 
the (too, TO1/2) are now allowed at the 95% CL, as 
compared to the pre-LHC fit, indicates that the 



^The NUHMl analysis includes both the dedicated 
NUHMl sample, which is efficient for smaller values 
of (mo,mx/2)i the basic set of CMSSM points, 

which provide extra NUHMl sampling at larger values of 
(mo,mi/2). 

^The parameters of the pre-LHC best-fit CMSSM and 
NUHMl points given in Table [T] differ by up to 1 c from 
those given in [8]. These differences are caused primarily 
by changes in the data inputs. 



tension between (g — 2)^, favouring relatively low 
SUSY scales, and the direct search limits, favour- 
ing larger SUSY scales, has significantly reduced 
the sensitivity of the fits within the CMSSM and 
the NUHMl for constraining the SUSY param- 
eters. Since the values of the best-fit points 
are significantly higher (see Table [1]), and conse- 
quently the distribution towards higher val- 
ues of mo and TO1/2 is much flatter than in the 
pre-LHC case, the precise location of the 68% 
and 95% CL contours is less precisely determined 
than before. The narrower range of toq allowed 
in the NUHMl at the 68% CL, compared to the 
CMSSM, is due to the appearance of a 'pit' with 
a lower absolute value of that is attainable 
in this model thanks to its flexibility in recon- 
ciling the {g — 2)^ with the Xi LSP and other 
constraints by deviating from Higgs mass univer- 
sality. As in our previous analysis including 2010 
LHC and XENONIOO data [5], we flnd no distinct 
enhancement of the likelihood in the focus-point 
region at large mo. 

The absolute values of x^ at the best-fit points 
for the pre-LHC case and for the LHCi/fb data 
set using our standard implementations of the 
{g — 2)^1 and BR(6 — > 57) constraints are given 
in Table [1] Our updated analysis of the pre- 
LHC data set yields xVd.o.f. = 21.5/20(20.8/18) 
in the CMSSM and NUHMl, respectively, cor- 
responding to p-valucs of 37% and 29%. On 
the other hand, using the LHCi/fb data set, we 
find that the minimum values of x^ are signif- 
icantly larger than the numbers of effective de- 
grees of freedom in the fits, which are also shown 
in Table [B xVd.o.f. = 28.8/22 (27.3/21) for 
the CMSSM and NUHMl, respectively Q. Cor- 
respondingly, the best fits have significantly re- 
duced probability values, 15% in the CMSSM 
and - 16% in the NUHMl B If we combine the 



'''For technical reasons, the Tz constraint was not included 
in our previous fits, leading to changes of one unit in the 
numbers of effective degrees of freedom in the fits. 
^The p-values for the VCMSSM and mSUGRA are some- 
what smaller, as was found previously [8]: x'^/^-o.f. = 
31.2/23(32.5/23), respectively, corresponding to p-values 
of 12% and 9%. We do not discuss these models fur- 
ther, except to comment that the light-Higgs funnel re- 
gion found previously in these models is now excluded by 
ATLAS data [18], in particular. 
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Figure 1. The (mo,mi/2) planes in the CMSSM (left) and the NUHMl (right). In each plane, the best-fit 
point after incorporation of the LHCi/fi^ constraints is indicated by a filled green star, and the pre-LHC 
fit 0/ by an open star. The Ax^ ~ 2.30 and 5.99 contours, commonly interpreted as the boundaries of the 
68 and 95% CL regions, are indicated in red and blue, respectively, the solid lines including the LHCi/^ 
data and the dotted lines showing the pre-LHC fits. 



BR(6 sj) errors linearly instead of quadrati- 
cally, the values of decrease by 0.8(0.7) in the 
CMSSM and NUHMl, respectively, and the p- 
values increase modestly to 18% in both cases. 
On the other hand, if we drop the {g — 2)^ con- 
straint, we find xVd-o.f = 21.3/20(20.3/19), re- 
spectively, corresponding to p-values of 38% in 
both cases. Thus, the qualities of our best fits 
are not very sensitive to the treatment of BR(6 
57), but are much more sensitive to the inclusion 
of (g — 2)^, as we discuss later in more detail. 

The degrees of non- universality, r = rn^^/rn^, 
for the best-fit NUHMl points in Table [1] are as 
follows: r ^ -57 (pre-LHC), r = -54 (LHCi/fb), 
r = —54 (linear BR(& — > 57) error combina- 
tion), r = —0.39 (dropping [g — 2)^ constraint), 
r = —0.39 (including both variants). Since r is 
quite poorly constrained, we do not quote its un- 
certainties. 

We note that the best-fit values of Mh are sig- 
nificantly higher in the CMSSM and NUHMl 
fits dropping {g — 2)^, with their large values 
of (too, TO1/2), than in the fits that include the 
{g — 2)^ constraint. An LHC measurement of 
Mh could provide a diagnostic discriminating be- 
tween models with light and heavy spectra of 



third-generation squarks, and help forecast the 
likelihood of discovering these sparticles in future 
LHC runs. Within the CMSSM and NUHMl, 
measuring Mh could advise us whether to take 
{g — 2)^ at face value or, conversely, hint towards 
extensions of those models. 

Table [2] displays the contributions to the to- 
tal of each of the observables at the best- 
fit points in the CMSSM and NUHMl, reveal- 
ing where this tension originates. We note, as 
already mentioned above, that there is an impor- 
tant contribution to the function coming from 
the (g— 2)^ constraint, as well as the LHCi/fb con- 
straints. We return later to a more complete dis- 
cussion of the tension between these constraints, 
and also of the treatment of the BR(fe — !> 57) con- 
straint. 

Table [5] also displays a similar break- 
down evaluated assuming the central values of 
muAa^il^{Mz) and Mz for a CMSSM point 
at very large (too,TOi/2) = (15,15) TeV, = 
100 GeV and tan/3 = 10, near the limit in which 
its predictions coincide with those of the SM for 
the same value of Mh- This is very similar to 
the global best-fit value of Mh in the SM ob- 
tained incorporating the limits from the direct 
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126.6l?;9 


Both 


19.5/19 


43% 


2020^ 


1410t 


25801 


48l 


126.61?;^ 



Table 1 



Comparison of the best-fit points found in the CMSSM and NUHMl pre-LHC (including the upper limit 
on BR(_Bs — > /i+/i^) available then), and with the LHCi/fo data set (also including the XENONIOO 
constraint) using the standard implementations of the [g — 2)^and BR(b — > sj) constraints, followed by 
variants first adding linearly the theoretical and experimental errors in BR(b —> S7J and then dropping [g— 
2)^, and finally combining both variants. The errors for fits dropping (_g — 2)^ are large and asymmetric, 
and are not indicated. The predictions for Mh do not include the constraint from the direct LEP Higgs 
search, and have an estimated theoretical error o/±1.5 GeV. 



Higgs searches at LEP, the Tevatron and the 
LHC [59]. Using the SM limit of the CMSSM 
within the MasterCode framework ensures that 
this evaluation of in the "SM" can be com- 
pared directly with those at the best-fit points 
in the CMSSM and NUHMl. In the "SM" case 
we discard the constraints imposed by the cos- 
mological dark matter density and XENONIOO, 
since there is no way to explain dark matter with 
the SM. We also discard the LHC missing-energy 
and H /A constraints, but all other constraints are 
kept. Consequently, we do list the contribution to 
from {g — 2)^. If this is included (omitted), the 
global for the "SM" is 32.7 (21.5). The num- 
ber of degrees of freedom for the "SM" is con- 
sequently 23 (22) and the p- value is 9% (49%). 
We observe that the p-value for the CMSSM is 
rather larger than that for the "SM" if (g — 2)^ is 
included, though similar if (g — 2)^ is not included 



in the "SM" and MSSM analyses. 

We also note that one of the big contributors 
to the global functions for all the models is 
^fb(^), which contributes Ax^ ~ 7 to each of the 
fits, suppressing all their p-values. Prior to the 
LHC results, the CMSSM yielded a significant 
improvement to x^ , with the dominant contribu- 
tion coming from (g — 2)^ (^X^ = —10.8). Other 
contributing observables were Mw [l^x^ = —1.6) 
and A£(SLD) (Ax^ = -1.3), though Aa,{b) was 
somewhat worse in the CMSSM (Ax^ = 2.1). 
The post-LHC comparison with the "SM" is 
shown in Tabled Looking at the entries for the 
electroweak precision observables, the only sig- 
nificant change is now that for [g — 2)^{/S.x^ ~ 
—6.8), with all other observables showing changes 
|Ax^| < 1. Thus, when [g — 2)^is dropped as a 
constraint, the resulting x^ for the best-fit points 
in the CMSSM and "SM" are very similar. In 
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Observable 


Source 


Constraint 




Ax^ 


Ax^ 




Th./Ex. 




(CMSSM) 


(NUHMl) 


("SM") 


mt [GeV] 




173.2 ± 0.90 


0.05 


0.06 


- 




42, 


0.02749 ± 0.00010 


0.009 


0.004 


_ 


Mz [GeV] 


l44i 


91.1875 ± 0.0021 


2.7x10"'' 


0.26 


- 


Vz [GeV] 


[26] / [15 


2.4952 ± 0.0023 ± O.OOIsusy 


0.078 


0.047 


0.14 


(^had [nb] 


/ mi 


41.540 ±0.037 


2.50 


2.57 


2.54 


Ri 


[26] / [m 


20.767 ± 0.025 


1.05 


1.08 


1.08 




m 1 m 


0.01714 ±0.00095 


0.72 


0.69 


0.81 




m 1 m 


0.1465 ± 0.0032 


0.11 


0.13 


0.07 


-Rb 


[26] / [S] 


0.21629 ± 0.00066 


0.26 


0.29 


0.27 




m 1 m 


0.1721 ± 0.0030 


0.002 


0.002 


0.002 


^fb(fc) 


[26] / [44] 


0.0992 ± 0.0016 


7.17 


7.37 


6.63 






n 0707 -\- no'^f^ 

U.U/U/ in U.UUoO 


U.oD 





so 




|26l / |44| 


0.923 ± 0.020 


0.36 


0.36 


0.35 


Ac 


[26] / [44] 


0.670 ± 0.027 


0.005 


0.005 


0.005 


^^(SLD) 


[26] / [44] 


0.1513 ± 0.0021 


3.16 


3.03 


3.51 


sin^ C(Qfb) 


[26] / [44] 


0.2324 ± 0.0012 


0.63 


0.64 


0.59 


M»- [GeV] 


EH] / Ha 


80.399 ± 0.023 ± O.OIOsusy 


1.77 


1.39 


2.08 




53; / WM, 


(30.2 ± 8.8 ± 2.0susy) x 10"'" 


4.35 


1.82 


11.19 (N/A) 


A/ft [GeV] 


[28] / [^ 56; 


> 114.4[±1.5susy] 


0.0 


0.0 


0.0 




HS] / [S] 


1.117 ±0.076exp 


1.83 


1.09 


0.94 




±0.082sM ± 0.050SUSY 








BR(B, ^ M'^M") 


[29] / m 


CMS & LHCb 


0.04 


0.44 


0.01 




m 1 m 


1.43 ± 0.43EXP+TH 


1.43 


1.59 


1.00 


BR(Bd 


m 1 m 


< 4.6[±0.01susy] X lO""* 


0.0 


0.0 


0.0 




mi m 


0.99 ± 0.32 


0.02 


< 0.01 


<0.01 




m 1 m 


1.008 ± 0.014exp+th 


0.39 


0.42 


0.33 




mi m 


< 4.5 


0.0 


0.0 


0.0 


^^^|XP/SM 


[49] / [5l][52] 


0.97 ± O.OIexp ± 0.27sM 


0.02 


0.02 


0.01 




[29] / [46][5l][52] 


1.00±0.01exp ±0.13sM 


■CO.Ol 


0.33 


<0.01 




W / fsT.'SJ 


1.08±0.14exp+th 


0.27 


0.37 


0.33 




[31] / m 


0.1120 ± 0.0056 ± 0.012susY 


8.4x10"* 


0.1 


N/A 




m 


(m^OjCTp') plane 


0.13 


0.13 


N/A 


jets + 


D1[2Q] 


(mo, mi/2) plane 


1.55 


2.20 


N/A 


H/A,H^ 


_;2t; 


(A'fA,tan/3) plane 


0.0 


0.0 


N/A 


Total xVd.o.f. 
p- values 


All 


All 


28.8/22 
15% 


27.3/21 
16% 


32.7/23 (21.5/22) 
9% (49%) 



Tabic 2 



List oj experimental constraints used in this work, including experimental and (where applicable) theoreti- 
cal errors: supersymmetric theory uncertainties in the interpretations of one-sided experimental limits are 
indicated by [...]. Also shown are the contributions that these constraints make to the total functions 
at the best-fit points in the CMSSM and NUHMl, respectively, and (for comparison) in the SM limit of 
the CMSSM (called "SM") including (excluding) (g — 2)^. The total values ofx^, the numbers of degrees 
of freedom and the p-values at these points are shown in the two bottom rows. 
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the case of the post-LHC NUHMl, wc also see 
a large drop in relative to the "SM" due to 
{g — 2)^(Ax^ = —9.4) and again all others give 
|Ax^| < 1. We also note that in both the CMSSM 
and NUHMl the best fits receive significant con- 
tributions from the LHC + jets searches. 

To illustrate further the impact of LHCi/fb ex- 
perimental constraints relative to pre-LHC pre- 
ferred regions, we display in Fig. [2] colour-coded 
contours of approximatqjp- values from our global 
fits for the CMSSM and NUHMl. Care is taken 
to count the effective number of degrees of free- 
dom at each point, considering all constraints 
that contribute non-trivially to the functions. 
Thus, for example, we drop the contribution 
of the LHCi/fb missing-energy constraints where 
they contribute Ax^ < 0.1, causing the visible 
changes in shading along drooping diagonal lines 
in both panels of Fig. [2l (This cut is applied only 
to the LHCi/fb missing-energy constraint.) Sub- 
stantial non-zero p- values are observed to extend 
to high mo and mi/2i in both pre- (upper panels) 
and post-LHC i/fb (lower panels), and both the 
CMSSM (left panels) and NUHMl (right pan- 
els) models. As also seen earlier, the primary 
effect of the LHCi/fb searches for jets + is 
most evident for mi/2, preferring higher values 
than that predicted by the pre-LHC global fits. 
At even higher (mo, mi/2), beyond the droop- 
ing diagonal line, slight increases in approximate 
p- values appear when comparing the pre-LHC re- 
suhs with the post-LHCi/fb and post XENONIOO 
results. This is due partly to the experimental 
constraint on BR(_Bs — > ^+/i~), which is nearing 
the SM prediction. Regions of the CMSSM and 
NUHMl parameter spaces approaching the high- 
mass decoupling limit receive a non-zero contribu- 
tion from BR{Bs — > A''''M~) in the post-LHCi/fb 
era, resulting in Ax^ < 1, which actaually im- 
proves the overall per effective degree of free- 
dom. Additionally, the XENONIOO constraint 
slightly prefers high mass scales, so as to accomo- 



^Strictly speaking, transforming a value to a p-valuc, 
using a specified number of degrees of freedom, is valid 
for Gaussian-behaved constraints. Because some of the 
experimental limits are one-sided and modelled in a non- 
Gaussian manner as previously described, the p- values re- 
ported here can therefore only be considered approximate. 



date the small "excess" in events, also resulting 
in slightly better values of pcr effective degree 
of freedom. 

Complementing the comparison of the p- values 
of the "SM", the CMSSM and the NUHMl, we 
now use the standard F-test to test the utility 
of adding one or several parameters to a model 
fit of data. Given a set of data comprising N 
observables and a model using m parameters, one 
may compute x^("^) for N — m degrees of freedom 
as done above. In general, adding r parameters 
produces a reduced value of x^(m + r), and the 
difference between these two x^ distributions is 
itself a x^ distribution for r degrees of freedom. 
The F-statistic is defined by 

Fx - + ^ > 0. (1) 

^ X^{m + r)/{N -m-r) ^' 

The probability that introducing the r new pa- 
rameters are warranted is found by integrating 
the F-distribution, pp{f,r,N — m ~ r), from 
/ ~ to 00. We use the F-test to illustrate 
the relative preference for various models. 

In our case, for the "SM" we have x^ ~ 
32.7(21.5) for 23 (22) degrees of freedom ff (.9-2)^ 
is included in (omitted from) the fit. Using the 
CMSSM value of x^ = 28.8 for 22 degrees of 
freedom, we find = 2.98, and the probabil- 
ity that switching to the CMSSM is warranted is 
Pf = 90%. Correspondingly, using the NUHMl 
value of — 27.3 for 21 degrees of freedom, we 
find F^ ~ 4.15, and the probability that switch- 
ing to the NUHMl is warranted is pp — 97%. We 
can also compare the improvement in x^ gained 
by moving from the CMSSM to the NUHMl. In 
this case the probability that the extra parame- 
ter needed to define the NUHMl model is pre- 
ferred over the CMSSM case is 71%. Fig. [3] uses 
shading to display values oipp in the (mo,mi/2) 
planes of the CMSSM and NUHMl. On the 
basis of these plots, we find that present data 
may warrant switching from the "SM" to the 
CMSSM or NUHMl for values of mi/2 out to 
~ 1500 GeV0. Beyond this range of mi/2, the 

-'''We note, however, that these results may be too 
favourable to the CMSSM or NUHMl, since they do not 
include the impacts of the many lower-sensitivity con- 
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m„ [GeV/c^] m„ [GeV/c^] 




Figure 2. The (mo,mi/2) planes in the CMSSM (left) and the NUHMl (right), for the pre-LHC data set 
(upper) and LHC'i/[\, data set (lower). In each plane, different regions are colour-coded according to the 
p-values found in our global fits. We note that in the LHCi/ny analysis the regions with p > 0.05 extend 
up to mi 12 2000 GcV in each case. 



motivations for these models would be signifi- 
cantly reduced. We also note that the F-test indi- 
cates that there would be no advantage in switch- 
ing to the CMSSM or the NUHMl if {g-2)^ were 
to be dropped. 

An important part of the motivation for low- 
scale supersymmetry is to alleviate the fine- 
tuning of the Higgs mass parameter in the Stan- 
dard Model. However, the problem of fine-tuning 
returns if the supcrsymmctric mass scales become 



straints from CMS and ATLAS. This problem could be 
avoided if the Collaborations publish official combinations 
of the sensitivities of their searches. 



large |60] . The required amount of fine-tuning is 
increased significantly in our LHCi/fb fits com- 
pared to our pre-LHC fits, principally because 
of the increases in the best-fit values of mg and 
7711/2- Specifically, in the CMSSM our best pre- 
LHC fit required fine-tuning by a factor ~ 100, 
whereas our best LHCi/fb fit requires fine-tuning 
by a factor ^ 300. The corresponding numbers 
for the NUHMl arc - 250 pre-LHC and - 600 
with the LHCi/fb data. 

Uncertainties in the analysis 

In assessing the compatibility of the CMSSM 
and the NUHMl with the experimental data it 
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Figure 3. The (fTiO:"^i/2) planes in the C'MSSM (left) and the NUHMl (right). In each plane, different 
regions are colour-coded according to the values of pp found by applying the F-test to our global fits. We 
note that the regions with pp > 0.8 extend up to mi^2 ~ 1500 GcV in each case. 



is important also to examine carefully the most 
important systematic uncertainties in the con- 
straints that drive the fit in this global likeli- 
hood analysis. We saw in the Table [2] that the 
most important contributions to the global 
functions at the best-fit points in these mod- 
els originate from the {g — 2)^ and LHCi/oj 
constraints. Exploring this further, Fig. |4] dis- 
plays the CMSSM and NUHMl (mo, TO1/2) planes 
again, exhibiting the contribution to from 
the {g — 2)^ constraint, evaluated for the model 
parameter sets that minimize the total 
each point in the plane. Their shapes are dif- 
ferent in the CMSSM and NUHMl, reflecting 
the existence of the previously-mentioned 'pit' 
in the latter model where x^ niay be reduced 
by some judicious choice of the degree of Higgs 
non-universality. Prima facie, there is tension 
between the (g — 2)^ constraint, which prefers 
small values of (?Tio,mi/2), and the LHCi/tb con- 
straints, which prefer larger values of (mg, mi/2). 
This tension is alleviated for larger values of 
tan/3, which is why post-2010-LHC [5] and post- 
LHCi/fb fits have favoured larger values of tan/3 
than pre-LHC fits [5] , albeit with large uncertain- 
ties. 

Fig. [5] again displays the (7710,7711/2) planes in 
the CMSSM (left) and the NUHMl (right), this 



time showing as solid lines the 68% and 95% 
CL contours obtained by dropping the {g — 2)^ 
constraint, the contours obtained applying the 
[g — 2)p constraint as in Fig. [1] being shown here 
for comparison as dotted lines. We see that, in 
the absence of the {g — 2)^ constraint, the outer 
parts of the 68% CL contours are expanded out- 
wards, close to the 95% CL contours that are 
themselves close to the boundary set by the ri^,o 
constraint 0- Within the overall range allowed 
by this constraint, the most important constraint 
is that provided by the LHC data. We note that 
the global likelihood functions in the CMSSM and 
NUHMl are very fiat, and that the best-fit points 
found dropping the {g — 2)^ constraints are cor- 
respondingly quite uncertain. However, it is clear 
that the amounts of fine-tuning at (g — 2)^-less 
best-fit points are much higher than if (g — 2)^ is 
included. 

We also comment on the treatment of BR(5 — > 
57). where different points of view have been 
taken concerning the combination of the theo- 



^ This constraint is not sacrosanct, but could be relaxed 
by postulating some amount of R-violation, or some other 
source of dark matter, or by modifying the expansion his- 
tory of the Universe, e.g., by altering the expansion rate 
during freeze-out, or by postulating some subsequent in- 
jection of entropy. 
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Figure 4. The (mo,TOi/2) planes in the CMSSM (left) and the NUHMl (right), with shading displaying 
the contribution to the global ^ function from (g — 2)^ as calculated using low-energy e^e~ data to 
evaluate the SM contribution. These contributions are evaluated for the parameter sets that minimize 
at each point in the planes. 




Figure 5. The (jtiq, mi/2) planes in the CMSSM (left) and the NUHMl (right). Here we show as solid 
lines the 68% and 95% CL contours obtained by dropping the {g — 2)^ constraint. The contours obtained 
applying the {g — 2)^ constraint as in Fig. Q] are shown here for comparison as dotted lines. 



retical and experimental errors [il- The two un- 



^^As in our previous analyses, we do not take into ac- 
count constraints from exclusive b — >■ 57 transitions. In 
particular, we do not impose any constraint on the SUSY 
parameter space from the isospin asymmetry in S — >■ K*j, 
as included for instance in the Superlso package [30) . A 
conservative treatment of non-factorizable contributions 
to this observable suggests a SM error exceeding ±0.05, 
i.e., a relative error exceeding 100% (see, e.g., I61| ). and 



certainties are of similar size, and the issue of 
how they are combined is more severe than for 
other observables. In our default implementa- 
tion of BR(6 — > s"f) we add the quoted errors in 



even larger errors are associated to the contributions of 
non-SM operators (see, e.g., 1621 ). These uncertainties ob- 
scure possible SUSY contributions within the ranges cur- 
rently of interest in the CMSSM and the NUHMl models. 
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quadrature. However, it might be argued that the 
theoretical error should be regarded as an overall 
range of uncertainty that cannot be treated as an 
effective Gaussian error to be added in quadrature 
to the experimental error. Therefore, we consider 
as an alternative implementation of BR(6 — >■ 57) 
the possibility of adding linearly the theoretical 
and experimental errors. As seen in Fig. [6l the 
region of the CMSSM {mo,mi/2) plane that is 
favoured at the 68% CL contracts significantly 
if the errors in BR(6 — > 57) are added linearly, 
whereas there is little effect in the NUHMl. 

This effect arises because the treatment of the 
BR(& — >■ 57) errors does not change the global 
function at large (toq 7^1/2)5 where its value ap- 
proaches that of the "SM" , since the experimen- 
tal measurement of BR(6 — >■ 57) is in good agree- 
ment with the SM prediction. On the other hand, 
adding the errors linearly relaxes the BR(6 — ?> 57) 
constraint at smaller (mo,TOi/2) in the CMSSM, 
reducing the tension with other observables and 
hence also the minimum of x^j seen in Ta- 
ble [TJ The net result is to enhance the rate of 
increase of for CMSSM parameters departing 
from the best fit, implying that the 68% CL is 
reached more quickly. On the other hand, this ef- 
fect is absent in the NUHMl because the freedom 
to adjust the degree of Higgs non-universality al- 
ready mitigates the tension of BR(& — > 57) with 
the other observables, leading to the 'pit' in the 
global function mentioned above. 

The {tan 13, 7711^2) planes in the CMSSM and 
NUHMl 

Fig. [7] displays the (tan/3, mi/2) planes in the 
CMSSM and NUHMl, exhibiting clearly the 
movement of the best-fit points and 68% and 
95% CL contours to larger tan /3 that is driven by 
the tension between — and the LHC push to 
larger mi/2- Comparing the dotted and solid con- 
tours, we see that the LHC Ipr constraints force 
the new best-fit points into what were previously 
the 'tails' of the 95% CL regions at large mi 12 
and hence tan /?. However, it is clear that the 
range of tan (3 allowed at the 68% CL is still very 
broad, extending from < 20 to > 50 in both the 
CMSSM and the NUHMl. On the other hand, 
any future substantial increase in the LHC lower 



limit on mi/2 would push tan/3 in both models 
into a narrower range ^ 50, where it encounters 
pressure from BR(_Bs — > M^A* ) discussed be- 
low. 

The (M^,tan/3) planes in the CMSSM and 
NUHMl 

We now turn to the (M^,tan/3) planes of the 
CMSSM and NUHMl, shown in Fig. El which are 
affected directly by the new CMS constraints on 
the heavy MSSM Higgs bosons H/A, H^ [I3[2T1 
SO], and by the CMS [22] and LHCb [23] con- 
straints on BR(_Bs As already dis- 
cussed, we include the heavy Higgs constraints 
via the most sensitive CMS search for H/A — > 
T+T~ |21| . In evaluating the constraint on the 
(M^,tan/3) planes of the CMSSM and NUHMl 
imposed by the new measurements of BR(_Bs — >■ 
/x^/x"), we use the official combination of the 
recent upper limits on BR(i?s — fi^fi^) from 
the CMS and LHCb Collaborations [IT], which 
yields BR{Bs ^ M+Z^") < 1.08 x 10"*^ at the 
95% CL. Our implementation actually includes 
the full likelihood function arising from this com- 
bination. The measurement of BR(i?s — >■ /Lt+/x~) 
by the CDF Collaboration [53] is in some ten- 
sion with the CMS/LHCb combination, but only 
at the Ax^ ~ 1 level. However, since there is 
so far no official combination of the CDF result 
with those of CMS and LHCb, we limit ourselves 
to discussing later its compatibility with the pre- 
dictions for BK{Bs — > of our global fits. 
We see in Fig. [8] that the general effect of the 
LHCi/fb data is to push the preferred range of 
Ma to larger values, as well as pushing tan /3 to- 
wards larger values. 

The impacts of the H/A and BR{Bs ) 
constraints are less important in the CMSSM 
than in the NUHMl, so we discuss the latter in 
more detail. Fig. [9] shows four versions of the 
(Af^,tan/3) plane in the NUHMl, with all the 
LHCi/fb constraints applied (upper left, equiva- 
lent to the right panel of Fig. [8]), dropping the 
H/A constraint but keeping BR(i?s — > fi^fi~) 
(upper right), dropping BR(i?s — !> but 
keeping the H/A constraint (lower left), and 
dropping both constraints (lower right). Compar- 
ing the two upper panels, we see that the H/A 
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m„ [GeV/c^] m„ [GeV/c^] 

Figure 6. The {111,0,7111/2) planes in the CMSSM (left) and the NUHMl (right). Here we show as solid 
lines the 68% and 95% CL contours obtained by adding linearly the experimental and theoretical errors 
in BR(b — > s"f). The contours obtained by combining them quadratically as in Fig. Q] are shown here for 
comparison as dotted lines. 




Figure 7. The (tan/3, mi/2) planes in the CMSSM (left) and the NUHMl (right). In each plane, the 
best-fit point after incorporation of the LHCi /fb constraints is indicated by a filled green star, and the pre- 
LHC fit by an open star. The Ax^ ~ 2.30 and 5.99 contours, commonly interpreted as the boundaries 
of the 68% and 95% CL regions, are indicated in red and blue, respectively, the solid lines including the 
LHCi/pa data and the dotted lines showing the pre- LHC fits. 



constraint is relevant for Ma 450 GcV, and 
that applying it impacts the low-M^ sides of the 
68% and 95% CL contours, whereas the best-fit 
point is unaffected. Comparing the upper and 
lower panels, we see that the BR(i33 — > /-t^yU^) 
constraint is relevant for low Ma values. In par- 



ticular, the 68% CL contours extend to slightly 
lower tan/3 values and the best-fit points (green 
stars) move to significantly lower tan /3 when 
BR(i3s — > /i^/-t^) is included, as seen in Fig. [3] 
Both the H/A and ^Ji{Bg — > ^+/Lt~) constraints 
have the potential for more significant impacts in 
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Figure 8. The {Ma, tan /?) planes in the CMSSM (left) and the NUHMl ( right). In each plane, the best-fit 
point after incorporation of the LHCi /fb constraints is indicated by a filled green star, and the pre-LHC 
fit 0/ by an open star. The Ax^ = 2.30 and 5.99 contours, commonly interpreted as the boundaries of the 
68 and 95% CL regions, are indicated in red and blue, respectively, the solid lines including the LHCi /fa 
data and the dotted lines showing the pre-LHC fits. 



the future. 

Predictions for rUg 

In Fig. [10] we show the one-dimensional 
functions predicted by global fits for rUg in the 
CMSSM (left) and the NUHMl (right). The solid 
lines are based on our global fits including the 
LHCi/fb constraints, whereas the dotted lines are 
based on our previous global fits based on the pre- 
LHC constraints [5] lij. Wc see that the best-fit 
estimates of irig have increased substantially to 
^ 1600 GeV as a result of the LHCi/fb data, but 
we also see that there is considerable uncertainty 
in this estimate, with nig > 2500 GeV subject to 
a penalty Ax^ ~ 2 only. 



l^The 'stalactites' at nig ~ 400 GeV in the pre-LHC 
fits [5] were due to the light-Higgs funnel that has now 
been excluded by the LHCi/fb data. Likewise, the 'sta- 
lactite' in the CMSSM LHCi/fb curve also originates in 
the light-Higgs funnel region, and comes from points with 
large mo > 3 TeV - which is why they are not seen in Fig.[l] 
- and large Ag. These points might be excluded by the 
ATLAS 1/fb 0-lepton search, whose published (mo,mi/2) 
exclusion for tan /3 = 10 and Aq = extends only to 
mo = 3 TeV. 



Predictions for BR(i?s — !> fi^fi^) 

In Fig. [TT] we show the one-dimensional 
functions predicted by our global fits for 
BR(B, ^ /-i+M") in the CMSSM (left) and the 
NUHMl (right). The solid lines are based on the 
official combination of the CMS and LHCb con- 
straints on this decay [H], whereas the dashed 
lines show results using an unofficial combina- 
tion of these constraints with the CDF measure- 
ment [21] , and the dotted lines represent pre-LHC 
predictions [3]. We see that the best- fit esti- 
mates of BR(i?s — >• fj,~^n~) are somewhat above 
the SM value, as a result of the push towards 
larger tan /3 required to accommodate the LHC 
data while reconciling them with (g — 2)^. In 
both the CMSSM and the NUHMl, the estimates 
of BR(i3s — ^ fi'^ fj,~) are quite compatible with an 
unofficial combined fit to CDF, CMS and LHCb 
data, where the main effect is a reduction of x^ in 
a somewhat broader range of BR{Bs — > fi'^^~). 

Predictions for Mh 

In Fig. [12] we show the one-dimensional 
functions predicted by our global fits for Mh in 
the CMSSM (left) and the NUHMl (right). In 
this figure we do not include the direct limits from 
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Figure 9. The (AfA,tan^) planes in the NUHMl including both the H/A HJ/ and BR{Bs M+M") R?F 
constraints (upper left), dropping the H/A constraint but keeping the BR(i3s — ^ fj,'^ constraint (upper 
right), dropping BR(i?s — fi'^i^i^) but keeping H/A (lower left), and dropping both constraints (lower 
right). 



LEP [SSlinn] or the Tevatron, so as to illustrate 
whether there is a conflict between these limits 
and the predictions of supersymmetric models. 
For each model we display the new likelihood 
functions corresponding to the LHCi/tb data set, 
indicating the theoretical uncertainty in the cal- 
culation of Mh of ^ 1.5 GcV by red bands. We 
also show, as dashed lines without red bands, our 
previous predictions based on the pre-LHC re- 
sults (also discarding the LEP constraint). We 
see that the LHC data improve the consistency of 
the model predictions with the LEP exclusion, re- 
moving whatever tension existed previously. We 



cannot resist pointing out that the best-fit value 
for ]\Ih found recently in a SM fit including LEP, 
Tevatron and LHC exclusions as well as precision 
electroweak data ~ 120 GeV [59], and that this is 
also the value of the SM Higgs mass that is most 
compatible with the ongoing LHC searches [55] . 

Predictions for Ma 

In Fig. [13] we show the onc-dimcnsional 
functions predicted by our global fits for Ma in 
the CMSSM (left) and the NUHMl (right). We 
see that the best-fit values of Ma have increased 
in both models, by - 350 GeV and ~ 250 GeV, 
respectively. 
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Figure 10. The one- dimensional functions for rrig in the CMSSM (left) and the NUHMl (right). The 
solid lines are for fits including the LHCi/^ data, and the dotted lines are for fits based on the pre-LHC 
data 0/. 
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Figure 11. The one- dimensional functions for IiB.{Bs in the CMSSM (left) and the NUHMl 

(right). The solid lines are for fits including the official combination of the results from the CMS and 
LHCb Collaborations 1^1^ , the dashed lines are for fits using an unofficial combination of these results 
with the CDF result \24^ , and the dotted lines represent pre-LHC predictions J^. 



Dark matter scattering cross sections 

In Fig. [H we show the 68% and 95% CL con- 
tours in the {m^a,ap^) planes for the CMSSM 
(left) and the NUHMl (right). The sohd lines 
are based on our global fits including the LHCi/fb 
constraints, whereas the dotted lines correspond 
to our previous fits using the pre-LHC con- 
straints. In both cases, we assume E^^r = 50 ± 



14 MeV [35] and we include with the LHCi/a 
data the XENONIOO constraint on [25]. We 
see that the LHCi/ft, data tend to push m-^o to 
larger values 0, and that these are correlated 



^ We recall the sensitivity of predictions for (7^ to the 
uncertainty in S^jv El- 

'^^The slivers of points at m-o ~ 60 GeV originate in the 
^1 

light-Higgs funnel region with large mo > 3 TcV men- 
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Figure 12. The one- dimensional functions for Mh in the CMSSM (left) and the NUHMl (right). The 
solid lines are for fits including all the available data, except the LEP l55\!tE^ constraints, with a red hand 
indicating the estimated theoretical uncertainty in the calculation of Mh of ^ 1.5 GcV. The pre-LHC 
results fSjl are shown as dotted lines. The yellow shading in the left panel shows the LEP exclusion of 
an SM Higgs boson, which applies also to the lightest CMSSM Higgs boson \5'7\\58^ . The lighter yellow 
shading in the right panel reflects the fact that this mass range is not completely excluded in the NUHMl 
due to a possible suppression of the ZZh coupling. The beige shading in both panels indicates values of 
Mh inaccessible in the super symmetric models studied with GUT-scale unification. 
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Figure 13. The one- dimensional functions for Ma in the CMSSM (left) and the NUHMl (right). The 
solid lines are for fits including the LHCx/q^ data, and the dotted lines are for fits based on the pre-LHC 
data l^. 

with lower values of though with best- fit val- 

tioned earlier, which might be excluded by the ATLAS ues Still ~ 10""*^ Cm^. We do not present here pre- 

1/fb 0-lcpton search. 
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dictions for spin-dependent scattering or signa- 
tures of astrophysical dark matter annihilations, 
which are further removed from the prospective 
experimental sensitivities in the near future. 

Sparticle thresholds in e^e~ annihilation 

In view of the interest in building an e^e~ 
collider as the next major project at the en- 
ergy frontier, we now analyze the implications 
of the LHCi/fb and XENONIOO data for expec- 
tations for sparticle production in e^e~ anni- 
hilation within the CMSSM and NUHMl. In 
this respect it has to be kept in mind that the 
LHC searches are mainly sensitive to the produc- 
tion of coloured particles, whereas lepton colliders 
will have a high sensitivity in particular for the 
production of colour-neutral states, such as slep- 
tons, charginos and neutralinos, as well as yield- 
ing high-precision measurements that will pro- 
vide indirect sensitivity to quantum effects of new 
states. Anything inferred from the coloured sec- 
tor concerning the uncoloured sector depends on 
the underlying model assumptions, and in partic- 
ular on assumptions about the possible universal- 
ity of soft supersymmetry breaking at the GUT 
scale. Non- universal models, e.g., low-energy su- 
persymmetric models, or models with different 
GUT assumptions, could present very different 
possibilities. 

Fig. [TS] compares the likelihood functions 
for various thresholds in the CMSSM (upper 
panel) and the NUHMl (lower panel), based 
on the global fits made using the LHCi/fb and 
XENONIOO constraints. The lowest thresholds 
are those for e+e" XiXiJ Tifi, eRCR and 
f^Rt^R (the latter is not shown, it is similar to 
that for J^eR). We see that, within the CMSSM 
and NUHMl, it now seems that these thresh- 
olds may well lie above 500 GeV, though in the 
CMSSM significant fractions of their likelihood 
functions still lie below 500 GeV. The thresholds 
for x5x2 ^'^'i Cijei -I- 6^6^ are expected to be 
somewhat higher, possibly a bit below 1 TeV. 
The preferred value for the threshold for xf 
lies at about 1700 GeV in both the CMSSM and 
NUHMl scenarios, that for the HA threshold 
lies above 1 TeV, and that for first- and second- 
generation squark-antisquark pair production lies 



beyond 2.5 TeV in both models. It should be 
kept in mind that these high thresholds are linked 
with the reduced p-value of the model. Fur- 
ther increases in the excluded regions would yield 
even higher thresholds, but would also make the 
CMSSM or NUHMl seem even less likely. 

4. Summary and Conclusions 

There is some disappointment in the air that 
the LHC has found no signs of supersymmetry 
in its first ~ 1/fb of data. However, it should 
be kept in mind that the searches performed at 
the LHC so far have essentially only been able 
to set limits on the production of the gluino and 
the squarks of the first two generations, and the 
resulting limits depend sensitively on the mass 
assumed for the lightest supersymmetric parti- 
cle [64]. On the other hand, the sensitivities of 
direct searches for stops and sbottoms and colour- 
neutral superpartners are very limited up to now. 
This situation will improve, as several times more 
data can be expected by the end of 2012, there 
is the prospect subsequently of an increase in the 
energy by a factor up to two, and the LHC is ex- 
pected eventually to accumulate orders of magni- 
tude more data. 

The initial optimistic prospects for SUSY 
searches at the LHC were largely driven by two 
indications that the supersymmetric mass scale 
might not be very high: [g — 2)^ and the need for 
dark matter that should not be overdense. Nei- 
ther of these indications has weakened recently. 
Indeed, the {g — 2)^ hint has even strengthened 
with the convergence of the previously discrepant 
SM calculations using low-energy e"^e^ and r de- 
cay data [HJinS]. However, as we have discussed 
in this paper, significant tension is now emerg- 
ing between the {g — 2)^ constraint and LHC 
data within the specific context of the CMSSM 
and NUHMl. A priori, in a general SUSY model 
there is not necessarily a tension between a heavy 
gluino and heavy squarks of the first two gener- 
ations on the one hand, as favoured by the LHC 
limits, and light colour-neutral states on the other 
hand, as favoured by {g — 2)^. 

The tension within the CMSSM and NUHMl 
can be reduced to some extent by adopting a 
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Figure 14. The 68% and 95% CL contours (red and 
NUHMl (right). The solid lines are for fits including 
the dotted lines include only the pre-LHC data JSjj. 

larger value of tan /3, but this may eventually 
lead to subsidiary tension with the LHC H/A 
constraints and the tightening experimental vise 
on BR(i?s — !■ In any case, it will be 

important to subject the {g — 2)^ constraint to 
closer scrutiny, and the upcoming Fermilab and 
J-PARC experiments on {g — 2)^ [66] are most 
welcome and timely in this regard. In parallel, re- 
finements of the experimental inputs for the pre- 
diction of (g — 2)^ from both low-energy e~^e~ 
and T decay data would also be welcome. It will 
be also necessary to subject the theoretical cal- 
culations within the SM and the corresponding 
estimates of the remaining theoretical uncertain- 
ties to further scrutiny. 

The dark matter upper limit on the sparticle 
mass scale remains unchanged, and is respon- 
sible for the disfavoured region above toi/2 ^ 
2500 GeV visible in our figures for the CMSSM 
and the NUHMl. On the other hand, the dark 
matter constraint on mo is not so strong, as also 
seen in the figures, extending well beyond the 
range displayed. Considering the impact of di- 
rect jets -I- l^T searches only, the regions of the 
CMSSM and NUHMl (mo,mi/2) planes in Fig.[2] 
with p-values significantly non-zero extend be- 
yond the likely reach even of the full-energy LHC 




blue, respectively) in the CMSSM (left) and the 
the XENONIOO JESI and LHC'i/it, data, whereas 



in its high-luminosity incarnation. A fortiori, the 
same is true for the regions of these planes allowed 
at the current 95% CL (Ax^ = 5.99 relative to the 
global minima, bounded by the blue contours in 
Fig. [T]). This is even more true of the full regions 
of the CMSSM and NUHMl (too,toi/2) planes 
that are allowed by the dark matter constraint. 

In light of this discussion, under what circum- 
stances could one conclude that the CMSSM or 
NUHMl is excluded? Currently, our best fits in 
both these models have p- values above 10%, com- 
parable to that of SM fits to precision electroweak 
data from LEP and SLD, and the F-test shows 
that both the CMSSM and NUHMl are war- 
ranted extensions of the SM, in the sense that in- 
troducing their parameters provides an improve- 
ment in that is valuable in both cases. More- 
over, it seems unlikely that the LHC will soon be 
able to explore all the region of the {mo,mi/2) 
planes in Fig. [2] where the models' p-values ex- 
ceed 5%, nor does the LHC seem likely soon to 
push Fy^ (see Fig. [3]) to uninterestingly low lev- 
els. This is not surprising, as in the high-mass 
limit the superpartners decouple and one is left 
essentially with the SM with a light Higgs. 

One way for the LHC to invalidate the mod- 
els studied here would be to discover an SM-like 
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Figure 15. The likelihood functions for various pair-production thresholds in e^e , as estimated in 
the CMSSM (upper panel) and the NUHMl (lower panel) after incorporating the XENONIOO JMI and 
LHCi/j\^ constraints. The likelihood function for the fiRfiB. threshold (not shown) is very similar to that 
for eRCR. 



Higgs boson weighing substantially more than the 
range 120 GeV predicted in Fig [121 A value of 
Mh 125 GeV or more would be in some tension 
with (g — 2)^, and perhaps hint towards mod- 
els beyond the CMSSM and NUHMl, whereas a 
value of ^ 130 GeV or more would cast severe 
doubt on most simple GUT-based models. As al- 
ready mentioned, range Mh ~ 120 to 130 GeV is 
precisely that currently favoured independently 
by precision electroweak data and by LEP, Teva- 
tron and LHC searches. If a Higgs- like signal were 
to be discovered in the lower part of this range, 
supcrsynimetry might not be far away, whereas 
if Mh is in the upper part of this range, indi- 
cating that at least the third-generation squarks 
could be heavy, one might for some time be in 
the frustrating situation of acquiring ever more 



circumstantial hints for supersymmetry, but with 
no direct evidence. On the other hand, if the LHC 
discovers that Mh > 130 GeV, the time might 
come to take another look at non-minimal super- 
symmetric models. 
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